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ABSTRACT 

We have assembled new Spitzer Space Telescope Infrared Array Camera observations of the 
mysterious binary star e Aurigae, along with archival far-ultraviolet to mid-infrared data, to form 
an unprecedented spectral energy distribution spanning three orders of magnitude in wavelength 
from 0.1 fim to 100 fim. The observed spectral energy distribution can be reproduced using a 
three component model consisting of a 2.2+q g M F type post-asymptotic giant branch star, and 
a 5.9±0.8 Mq B5±l type main sequence star that is surrounded by a geometrically thick, but 
partially transparent, disk of gas and dust. At the nominal HIPPARCOS parallax distance of 625 
pc, the model normalization yields a radius of 135 ± 5 Rq for the F star, consistent with published 
interferometric observations. The dusty disk is constrained to be viewed at an inclination of i > 87°, 
and has effective temperature of 550 ± 50 K with an outer radius of 3.8 AU and a thickness of 0.95 
AU. The dust content of the disk must be largely confined to grains larger than ~ 10 /im in order to 
produce the observed gray optical-infrared eclipses and the lack of broad dust emission features in 
the archival Spitzer mid-infrared spectra. The total mass of the disk, even considering a potential 
gaseous contribution in addition to the dust that produces the observed infrared excess, is <C 1 Mq. 
We discuss evolutionary scenarios for this system that could lead to the current status of the stellar 
components and suggests possibilities for its future evolution, as well as potential observational tests 
of our model. 
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1. INTRODUCTION 

The bright star e Aurigae (HD 31964) is a single-lined 
spectroscopic binary that is famous for its long orbital 
period (27.1 yr), which is punctuated by an almost 2 
yr long eclipse caus ed by an essentially invisible object 
(jCarroll et al.1 11991). The central problem posed by this 
system is that if the F star component, which dominates 
the light from the system over a wide range of wavelength 
and is the eclipsed object, is a massive supergiant (as its 
spectrum implies), then the invisible companion is sur- 
prisingly under-luminous for its mass. Exotic solutions 
for t his mass conund rum involving, for example, a black 
hole ()Cameronlll97lT ) are not viable because of the train 
of ever more complicated additional requirements that 
observational constraints impose on such a model. For 
example, the lack of significant X-ray emission from the 
system precludes a black hole unless there is no accretion 
from the disk, which is not possible unless there is yet 
another unseen body (a massive planet, perhaps?) that 
clears out the space around the black hole, and so on 
(see t he d iscussions inlBarsonv et al.lll986l iCarroll et al.l 
Il99ll and lWolk et aLll201ffl . 

By examining the optical spectra of e Aur near the 
end of its 1954-1956 eclipse, EacS dl959h was able to de- 
duce the electron density and develop the hypothesis of 
a Be star-like hot object at the center of a large disk of 
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occulting material pack! 119611 ). IWooU (|1973f ) reported 
pioneering infrared (IR) observations that revealed the 
presence of an excess consistent with the disk being a 
cloud of partially ionized gas, with a total projected area 
comparable to that of the F star. An estimate of elec- 
tron density from the IR excess was consistent with the 
optical-ultraviolet (UV) estimates of 10 1 1 cm" 3 . 

As IR detector technology advanced, [Backman et al.l 
(1984) reported ground-based IR photometry obtained 
during the 1982-1984 eclipse that demonstrated that 
the excess could be characterized a s a T = 500 ± 150 
K sou rce subtending 8 x 10~ 16 sr. Backma n fc Gillett] 
(|1985[ ) refined this result with Infrared Astronomy Ex- 
plorer (IRAS) satellite photometry during the eclipse, 
extending the wavelength coverage well into the thermal 
IR, and yielding a revised temperature of 475 ± 50 K and 
angular extent of 8.6 ± 1.0 X 10" 16 sr. iSticklandl (|1985f ) 
examined the same IRAS data and was led to conclude 
that the disk temperature could be better characterized 
as a 750 ± 100 K source with a projected area about 
six times that of the F star photosphere. Regardless of 
its exact characteristics, the transiting disk in the e Aur 
binary offers a valuable opportunity to study its longi- 
tudinal structure in ways not possible with circumstellar 
disks around single stars. 

We report here on the results from new mid-IR ob- 
servat ions of e Aur made with the Spitzer Space Tele- 
scope ([Werner et al.|[2004l) , which provide a more precise 
characterization of the occulting body, as well as a new 
look at archival data at shorter wavelengths, which bet- 
ter constrain the stellar components. The importance 
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of this exercise lies in the fact that the spectral energy 
distribution (SED) now can be much more precisely de- 
fined, thanks to the availability of new and recalibrated 
data spanning the far-UV to the mid-IR. These results 
strongly imply that the putative F supergiant star in 
e Aur is more likely a lower mass, unstable post-AGB 
object that previously transferred matter t o a B( e)-like 
star companion - as proposed bv lWebbink (11985) - re- 
sultin g in a complex "dark matter" disk ()Howell et all 
2008) that causes the eclipses. 

2. SPECTRAL ENERGY DISTRIBUTION - THE DATA 

To construct our SED for e Aur, which spans three 
orders of magnitude in wavelength, we combined data 
from new observations in the mid-IR from Spitzer, recent 
observations from the American Association of Variable 
Star Observers (AAVSO) in the optical, and archival 
ground- and space-based observations at other wave- 
lengths. The space-based data in particular come from 
the Midcourse Space Experiment (MSX), the Interna- 
tional Ultraviolet Explorer (IUE), the Hubble Space Tele- 
scope (HST), and the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE). The SED data are listed in Table [Hand 
plotted in Figure [TJ All of the observations were ob- 
tained outside of eclipse phases; the majority were ob- 
tained prior to the onset of the eclipse that began in 
late-2009, but well after the end of the prior eclipse in 
1984, except as noted in Table [TJ We describe our new 
Spitzer mid-IR observations in more detail below. 

2.1. Spitzer Infrared Array Camera 

The expected flux density of e Aur at 3-5 //m exceeds 
the nominal saturation limits of channels 1 (3.6 /jm) 
and 2 (4.5 jum) of the Infrared Array Camera (IRAC; 
iFazio et all 120041) on Spitzer by factors of several times. 
The rationale for observing e Aur with IRAC was two- 
fold: (1) as a test case for a potential monitoring cam- 
paign whose cadence could not be accomplished from the 
ground, and (2) to obtain measurements in the 3-5 /zm 
range whose flux calibration is homogenous with that of 
the other Spitzer data already obtained at longer wave- 
lengths (see Section 12.21 and l 2.3p . As it turned out, a 
third reason became apparent after the fact; namely, that 
existing older IR data in this wavelength range are not 
directly comparable to the more recent longer wavelength 
Spitzer data (see discussion in Section [3.4[1 . 

We utilized an allocation of six minutes^ of Director's 
Discretionary Time in late- April 2009, shortly before de- 
pletion of the Spitzer cryogen in mid-May 2009. We used 
a novel observing strategy that took advantage of both 
the extremely short exposure time available in the IRAC 
sub-array mode and the specific placement of the target 
centroid at the intersection of four pixels. The latter con- 
dition spreads the total illumination of the detector, as 
well as the brightest part of the point response function, 
over the maximum number of pixels. In sub-array mode, 
observations are obtained in units of 64 exposures of a 
32 x 32 pixel section of the full 256 x 256 pixel array, and 
are tiled onto the full array for downlink. This resulted 
in four dithered observations of e Aur in each channel, 

4 Three minutes of which was the standard "slew tax" applied 
to all Spitzer observations. 



with each dither comprised of 64 consecutive 0.02-sec ex- 
posures. 

We initially applied the IRAC array-location- 
dependence correction to the individual sub-array 
images, as described in the IRAC Data Handbook, and 
then performed aperture photometry on them (e.g., as 
described in IHoard et all 120091) using the IRAF0 task 
phot. We utilized a 10-pixel radius aperture (1 pixel « 
1.22"), with a 10-20 pixel background annulus, which is 
the configurations for which the corresponding aperture 
corrections are 1.00 in all IRAC channels. 

The high time resolution and photon-abundant nature 
of these observations revealed two interesting instrumen- 
tal artifacts that are normally below the detection thresh- 
old in Spitzer IRAC observations: a rapid "jitter" with 
amplitude of ~ 1% in the measured flux densities in 
each sequence of sixty-four 0.02-sec exposures, as well as 
overall trends with slopes of 1-2% in measured flux den- 
sity during each 64-frame sequence. The former effect 
is caused by sub-pixel mismatch between the measured 
and actual centroid of the point response function (PRF), 
propagated through the normal IRAC pixel-phase correc- 
tion, which accounts for the exact position of the target 
centroid within a pixel. It is likely caused by resolving 
the small discrete steps in the Spitzer tracking motion. 

The latter effect appears to be linked to the telescope 
settle time of 2 sec after a slew or dither offset, which 
is longer than the « 1.5-sec total length of one of the 
64-frame sub-array sequences. Both of these effects are 
presumably present during the first 2 seconds of any nor- 
mal, full-array IRAC exposure, for which the shortest 
available exposure time is 2 sec. Consequently, in a full- 
array observation of a fainter target, these effects are not 
time-resolved and are below the level of photon counting 
noise, so go unnoticed in the final, single photometric 
measurement per observation. In addition, we see the 
offsets of order several per cent between mean flux densi- 
ties at different dither positions (caused by pixel-to-pixel 
variations) that are also found in dithered full-array ob- 
servations, and are mitigated by averaging together the 
flux densities measured at the different dither positions. 

Additionally, in the channel 2 observations, we see a 
steep linear trend in the measured flux densities in each 
of the four dithered sequences of 64 exposures. This re- 
sults in an increase of «5-7%, with approximately equal 
slope from the start to the end of each dither sequence. 
This could possibly include a contribution from the so- 
called "ramp" effect noticed in so me exoplanet transit 
light curves obtained with IRAC (jDemi ng 2009), but, 
considering the very short total duration of these expo- 
sure sequences, is more likely due to build-up of latents 
on the pixels exposed to e Aur during each sub-array 
sequence. This particular artifact is not present in our 
channel 1 data. We note that a 2.5-hr sequence of deep 
IRAC exposures in the Oph star forming region ended 
about an hour before our e Aur observation, and im- 
mediately prior to our observations there was a short 
(20-min) series of observations in the Galactic plane. So, 
it is entirely possible that the charge traps in channel 
1 that would lead to latents were already full when our 

5 See http://ssc.spitzer.caltech.edu/irac/dh/ 

6 The Image Reduction and Analysis Facility is maintained and 
distributed by the National Optical Astronomy Observatory. 
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Figure 1. Observed (dereddened) SED of e Aur with a three component model. From short to long wavelengths, the photometric points 
are: U, B, V, R, I from the AAVSO (filled circles), J, H, K, L, M (filled circles), J, H, Ks from 2MASS (unfilled circles), IRAC from Spitzer 
(filled squares), ground-based L/, M, N, Q (unfilled squares), Bl, B2, A, C, D, E bands from MSX (unfilled diamonds), and MIPS from 
Spitzer (filled diamonds). Vertical error bars are the photometric uncertainties (which are dominated by the systematic uncertainty of the 
dereddening process for the dereddened data). The spectroscopic data are: FUSE (dark purple), ffi>T-GHRS (light purple), IUE (dark 
blue), ground-based optical (light blue and green), IRS (orange) and MIPS-SED (filled red squares) from Spitzer. See the text and Tablefl] 
for more information about the data. The model (solid line) is the sum of limb-darkened model F0 (post-AGB) and B5 V spectra (dotted 
lines), and a cool blackbody disk (dashed line). See the text and Table [T] for more information about the model. 



observation started. 

The final photometric measurements reported in Table 
[T]were obtained as follows: for channel 1 (IRAC-1), in or- 
der to mitigate against one saturated pixel in the e Aur 
PRF during the first dither sequence, we combined all 
256 sub-array images using the Spitzer data analysis soft- 
ware tool MOPEX (MOsaicker and Point source Extrac- 
tor )Q, and performed aperture photometry (as described 
above) on the mosaic. This results in a flux density that 
is 0.3 Jy fainter than the value obtained by averaging to- 
gether the aperture photometry results for all 256 indi- 
vidual sub-array exposures, including the one saturated 
pixel in 64 exposures. Qualitatively, this result is consis- 
tent with the removal of the saturated pixel in 64 of the 
images, resulting in a slightly smaller "true" flux density 
measurement. For channel 2 (IRAC-2), there are no sat- 
urated pixels in any of the individual sub-array images; 
however, in order to mitigate against the latent charge ef- 
fect described above, we averaged together only the first 
5 exposures in each dither sequence. In addition, we had 
to utilize a 5-pixel aperture (with 5-10 pixel background 
annulus) because the 10-pixel aperture was partially off 
the sub-array during two of the dither sequences. This 
required use of the corresponding aperture correction for 
IRAC-2 of 1.064. 

In all cases, as described in the IRAC Data Handbook, 

7 See http : //ssc . spitzer . caltech . edu/postbed/mopex . html 



we applied the pixel-phase correction (|Hora et~all[200l) 
to the photometry in both channels. However, we did 
not perform the color correction other than to utilize 
the isophotal effective channel wavelengths (see Table [TJ 
during our subsequent interpretation of the data, which 
accounts for almost all of the color correction. The re- 
maining effect of the color correction is folded into our 
IRAC system atic uncertainty bu dget. As described in 
more detail in iHoard et alj (|2009f ). the total uncertainty 
budget for IRAC photometry includes several systematic 
terms (3% for the absolute gain calibration, 1% for re- 
peatability, 3% for the absolute flux calibration of the 
IRAC calibration stars, and 1% for the color correction), 
as well as a scatter term evaluated from the r.m.s. scatter 
of the individual flux density measurements after remov- 
ing all long timescale trends. The uncertainties in the 
IRAC photometry values listed in Table Q] reflect these 
terms added in quadrature. 

2.2. Spitzer Infrared Spectrograph 

A Spitzer GO-2 observing program (20058; IStencel 
12007( 1 on e Aur obt ained observations w ith the Infrared 
Spectrograph (IRS; iHouck et all 12004( 1 and the Multi- 
band Imaging Photometer for Spitzer (MIPS; see be- 
low). For this work, we have utilized the current mature 
pipeline-(re)processed data for those observations, and 
re-extracted final data products. The two high resolu- 
tion IRS modules were used during two visits to e Aur, 



Table 1 

The Data 



Wavelength 


Band 


Source 


Date of Observation: 
UT JD 


Orbital 
Phase a 


Flux Density: 
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(■ r >) 

(6) 

8 

S 
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0.111-0.117 spectrum FUSE 010107 51917 0.646 0.00025-0.0015 0.15-0.3 Ake (20061 

0.117-0.146 spectrum HST-GHRS 960216 50130 0.466 0.002-0.02 0.4-1.2 T^ffeTXLambert (1999) 

0.150-0.198 spectrum IUE-SWP 850203, 850317 46099, 46141 0.058-0.062 0.005-0.8 0.6-20 Sheffer fe Lambert (1*9991 

0.185-0.335 spectrum IUE-LWP 861119-861123 46753-46757 0.124-0.125 0.5-48 10-130 Sheffer fc Lambert (1999) 

0.329-0.549 spectrum ground 820405 45065 0.953 25-215 150-630 Thompson et al. (1987) 

0.567-0.889 spectrum ground c. 1990-1992 c.47892-48987 0.25-0.33 240-330 660-520 Torres-Dodeen fc Weaver (1993) 

0.360 U AAVSO 031101-090701 52918-55013 0.748-0.959 59.1±1.8 306.0±37.6 this work 

0.44 B AAVSO 031101-090701 52918-55013 0.748-0.959 148.4±4.5 591.5±49.4 this work 

0.55 V AAVSO 031101-090701 52918-55013 0.748-0.959 230.4±6.9 663.9±59.3 this work 

0.70 R AAVSO 031101-090701 52918-55013 0.748-0.959 294.6±8.8 605.8±55.9 this work 

0.88 I AAVSO 031101-090701 52918-55013 0.748-0.959 378.4±11.4 606.8±67.9 this work 

1.20 J ground 970907-000418 50699-51652 0.523-0.620 297.1±10.3 386.9±53.4 Taranova fc Shenavrin (2001) 

1.235 J 2MASS 991108 51491 0.603 282.2±78.6 362.5±112.3 Skrutskie et al. (2006) ' 

1.62 H ground 970907-000418 50699-51652 0.523-0.620 254.8±14.0 296.8±47.5 Taranova fe Shenavrin (2001) 

1.662 H 2MASS 991108 51491 0.603 213.6±35.5 247.1±55.5 Skrutskie et al. (2006) ' 

2.159 Ks 2MASS 991108 51491 0.603 162.5±32.4 177.7±45.4 Skrutskie et al. (2006) 

2.20 K ground 970907-000418 50699-51652 0.523-0.620 177.1±6.1 193.1±31.7 Taranova fc"§henavrin (2001) 

3.50 L ground 970907-000418 50699-51652 0.523-0.620 94.5±3.3 98.0±16.9 Taranova fc Shenavrin' (2&tt) 

3.544 IRAC-1 Spitzer 090426 54948 0.953 66.2±3.0 68.6±12.0 this work 

3.8 L' ground 811113 44922 0.939 76.5±2.6 79.0±13.7 Backman et al. (1984) 

4.29 MSX-B1 MSX c.960424-970222 c.50198-50502 0.47-0.50 67.9±6.1 69.7±13.4 Egan et al.720*b3*' 

4.35 MSX-B2 MSX c.960424-970222 c.50198-50502 0.47-0.50 72.1±6.6 74.0±14.3 Egan et al. (2553) 

4.479 IRAC-2 Spitzer 090426 54948 0.953 52.9±2.4 54.2±9.6 this work 

4.8 M ground 800130-811210 44269-44949 0.873-0.942 56.8±2.9 58.0±10.4 Backman et al. (1984) 

5.00 M ground 970907-000418 50699-51652 0.523-0.620 56.3±3.4 57.4±10.5 Taranova fc"Si7enavrin (2001) 

8.28 MSX-A MSX c.960424-970222 c.50198-50502 0.47-0.50 24.4±1.0 ■■■ Egan et al. (2003) 

9.89-37.14 IRS-1, -3 Spitzer 051019, 060317 53663, 53812 0.823, 0.838 15.3-5.1 ■ • ■ this work 

10.1 N ground 800201-811217 44271-44956 0.873-0.942 14.5±0.5 ■•■ Backman et al. (1984) 

12.13 MSX-C MSX c.960424-970222 c.50198-50502 0.47-0.50 11.5±0.6 ■•■ Egan et al.720031 

14.65 MSX-D MSX c.960424-970222 c.50198-50502 0.47-0.50 9.6±0.6 ■•■ Egan et al. (2003) 

20.0 Q ground 811210-811217 44949-44956 0.942 3.8±0.2 ■•■ Backman et a!"Tl984) 

21.34 MSX-E MSX c.960424-970222 c.50198-50502 0.47-0.50 4.7±0.3 ■•■ Egan et al.725b3l 

23.675 MIPS-24 Spitzer 050925, 060223 53639, 53790 0.820, 0.836 3.7±0.2 ■ • ■ this work 

53.7-91.1 MIPS-SED Spitzer 050925, 060223 53639, 53790 0.820, 0.836 0.88-0.41 ■ • ■ this work 

71.44 MIPS-70 Spitzer 050925, 060223 53639, 53790 0.820, 0.836 0.5±0.07 ■ • ■ this work 



(7) 
(7) 
(8) 
(4,9) 
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(7) 

(4,9) 

(7) 

(1,10) 
(4,9) 



-: 

o 



(4,9) 

(1,11) 
(1,12,13) 

(1,11) 



Note. — (1) Flux densities are the averages of 2 measurements on the indicated dates. (2) IUE exposures SWP 25156 and 25470, and LWP 09565, 09547, 
and 09554. (3) Flux densities are the averages of 3 measurements in the indicated date range. (4) Note that the orbital phase refers to the previous cycle (i.e., 
before the 1982-1984 eclipse). (5) Exact date of spectrum is not given in the corresponding publication. (6) Flux densities are the averages of N data points in 
the indicated date range: Njj — 411, N& — 411, iVv — 411, — 37, Ni — 21. (7) Flux densities are the averages of 34 measurements in the indicated date range. 
(8) Spitzer AOR key 33903360, processing pipeline version S18.7.0. (9) Flux densities are the averages of N pre-eclipsc data points in the indicated date range: 
7V L / = 1, N M = 4, TVn = 4, N Q = 2. (10) Spitzer AOR keys 13848832 and 13849600, processing pipeline version S18.7.0. (11) Spitzer AOR keys 13849088 and 
13850112, processing pipeline version S16.1.0. (12) Spitzer AOR keys 13849856 and 13849344, processing pipeline version S16.1.0. (13) MIPS-SED flux densities 
have been normalized to the MIPS-70 photometric point by scaling the former by a factor of 0.84. 
a Using the orbital ephemeris JD oba = JD2445525 + 9890E (e.g., [Schmidtke et aT||1985l ; ICarroll et al.||1991|) . 
k See text for a discussion of the dereddening applied to the data. 

c When necessary, photometric data reported in magnitudes were converted to flux densities using appropriate filter-specific zero point values (e.g., see lCohen et all 
[2003T : IBessell fc Brettlfl988l ; ICampins eTal~lll985l ; IBessell|[T979l ; IBeckwith et al.|fl976l ; IMendoza Vl[l963T) . 

For this work, we have re-extracted all photometric and spectroscopic measurements from archival data, but list here the first original publication of the relevant 
data. 



e Aurigae from the Far-UV to the Mid-IR 



5 



and we utilized the Spitzer data analysis software tool 
SPICE (Spitzer IRS Custom Extraction j3 to extract the 
calibrated spectra from each visit for each of the two 
nod positions from the combined post-BCD data prod- 
ucts (which simply co-add the three individual cycle ex- 
posures at each nod position from each visit). Module 1 
spans 9.9-19.5 /im, while module 3 spans 18.8-37.1 /xm, 
at a resolving power of ~ 600 in both modules. Prior 
to extraction, the Spitzer data analysis software tool 
IRSCLEAJN0 was used to remove both permanent and 
rogue hot pixels from the combined post-BCD images. 
Because of the short exposure times used for individual 
frames (6- and 14-sec in modules 1 and 3, respectively), 
as well as the brightness of the target, no offset sky back- 
ground spectrum was subtracted; this is in keeping with 
the "best practices" recommendation for observations of 
this type from the Spitzer Observer's Manua{3- The four 
extracted spectra from each module (2 nods from each of 
2 visits) were then combined in a weighted average after 
rejecting all wavelength points flagged as unreliable by 
SPICE. 

2.3. Multiband Imaging Photometer for Spitzer 

The MIPS (|Rieke et al J 120041 ) data on e Aur were ob- 
tained in two observing modes. The first of these was the 
normal small field photometry mode at 24 and 70 /xm 
(MIPS-24 and MIPS-70, respectively). We performed 
aperture photometry on the filtered combined post-BCD 
ima ges following the procedures in the MIPS Data Hand- 
bool|3 operationally, this process is very similar to that 
used to obtain the IRAC photometry. For the MIPS-24 
images, we utilized a 35-arcsec aperture with a 40-50- 
arcsec background annulus, and applied the correspond- 
ing aperture correction of 1.08. For the MIPS-70 images, 
we utilized a 35-arcsec aperture with a 39-65-arcsec back- 
ground annulus, and applied the corresponding aperture 
correction of 1.22. As with the IRAC photometry, we did 
not apply the color correction because most of it is ac- 
counted for by utilizing the isophotal band wavelengths 
(see Table [lj. Instead, we folded the remaining small 
effect into our uncertainty budget. 

For MIPS-24 photometry, the total systematic uncer- 
tainty budget is the quadrature sum of 4% for the abso- 
lute calibration, 0.4% for repeatability, 3% for the color 
correction, and 5% for the aperture correction. These 
values reflect upper limits from version 3.3 of the MIPS 
Data Handbook. In addition, the total uncertainty in- 
cludes a scatter term of 0.006 Jy (0.2%) obtained for 
the target aperture in the uncertainty ( "munc" ) image 
provided as part of the post-BCD data products. For 
the MIPS-70 photometry, the corresponding terms in the 
uncertainty budget are: 7% for the absolute calibration, 
5% for repeatability, 3% for the color correction, 5% for 
the aperture correction, and a scatter term of 0.007 Jy 
(1.4%). 

The other MIPS observation utilized the MIPS-SED 
mode, which obtains a very low resolution (R «15— 25) 
spectrum between 52-97 /xm. We performed a weighted 
average of the post-BCD extracted spectrum from each of 



See http : //ssc . spitzer . caltech . edu/postbcd/spice . html 
See http: //ssc . spitzer . caltech. edu/postbcd/ irsclean.html 
5 See http://ssc.spitzer.caltech.edu/documents/SQM/ 
See http://ssc.spitzer.caltech.edu/mips/dh/index.html 



the two MIPS-SED visits to e Aur. The flux calibration 
uncertainty for MIPS-SED mode is ~ 20%, and we found 
it necessary to scale the data by a factor of 0.84 (i.e., 
a 16% reduction compared to the extracted values) in 
order to match the MIPS-SED spectrum with the better 
calibrated MIPS-70 photometric point. 

3. SPECTRAL ENERGY DISTRIBUTION - THE MODEL 

It has been fairly well establis hed (e.g., seelHuan j|1965l 
IKopallH97l1 and the review in IWebbink|[l985ft~ that the 
e Aur system consists of three primary components: an 
F star orbited by one or two B stars, with a large dusty 
disk surrounding the latter. It is the disk that eclipses 
the F star every 27 years. The exact details of these com- 
ponents, however, remain rather nebulous; for example, 
the mass of the F star ha s been proposed to b e either 
very high (10-20 M Q ; e. g.. iLissauer et all 1 1996ft or rela- 
tively low (1-4 Mn ; e.g. JLissauer et al.lll9 96: Sait o et all 
ll987tlTakeut3l986ft . The disk around the B star has been 
variously described as thick or thin, flat or twisted, edge- 
on or slightly inclined, opaque or semi-transparent, and 
"solid" or containing a central void. 

Fortunately, in the nearly three decades since its last 
eclipse, some important new information about e Aur has 
been obtained. This includes a HIPPARCOS trigono- 
metric parallax distance (with a nominal value of 625 
pc, which we have used i n our model calculations be- 
low; [Perrymair^Fjil][l993), various interferometric mea- 
surements of the angular size of t he F star (sa 2.1- 
2.3 mas from 0.45 /xm to K-band: iStencel et "ail 120081: 
iMozurkewich et al.ll200l iNordgren et al.ll2001ft 7 and the 
space-based far-UV and mid-IR observations shown in 
Figure [U which reveal details about the system compo- 
nents that are not the F star (which otherwise completely 
dominates the SED of e Aur in the range 0.2-4 /xm). Us- 
ing these data, we are now able to definitively constrain 
some of these heretofore uncertain system parameters. 
Our observed and model SEDs are shown in Figure Q] 
and the model parameters are listed in Table [T] We dis- 
cuss the model parameters, constraints, and components 
in more detail below. 

3.1. Dereddening of the Observed SED 

The observed flux densities were dereddened using the 
UV-o ptical-IR extinction law from lFitzpatrick fc Massal 
([2007ft . For A > 5 /xm, the reddening correction was 
negligible, so was not applied to the data. An in- 
terstellar reddening to e Aur, E(B - V) = +0.38, 
was used for A = 0.33-5 /xm. This value was fine- 
tuned from a starting value of E(B — V) = +0.35 in 
the r ange E(B — V) sa 0.3-0.4 found in the literature 
(e.g.. IMozurkewich et al.ll2003l:lTorres-Dodgen fc Weaverl 
[19931: iSticklandl 11985] lAkd 119851: Hack fc Selvellil Il979t 
iCastellil 119781 : iHobbs 1969), in order to best match the 
optical spectra and photometry in the SED to the model 
F star spectrum. We note that the apparent disagree- 
ment between the optical spectrum at the shortest wave- 
lengths, near 0.3 /xm, most likely reflects issues in accu- 
rately calibrating a ground-based UV spectrum, rather 
than an actual disagreement with the model, since the 
model matches the flux density of the IUE spectrum that 
ends at 0.3 /xm quite well. At wavelengths shortward of 
A = 0.33 /xm, additional dereddening was applied - see 
Section [331 for details. 
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Table 2 

The Model 



Component 



Parameter 



Value 



Reference 



System Adopted Distance, d (pc) 

Inclination, i (°) 
Orbital Separation, a (AU) 

F Star Spectral Type 

Temperature, Tp (K) 
log g 

Radius, Rp (Rq) 

Angular Diameter, D a (mas) 

Mass, M F (Mq) 

B Star Spectral Type 

Temperature, Tb (K) 
logfl 

Radius, i?B (Rq) 
Mass, M B (Mq) 

Disk Temperature, T^isk (K) 

Radius, i?disk (AU) 
Height, H disk (AU) 
Assumed Mass, M disk (Mq) 
Inferred Dust Grain Radius, r s 
Transmissivity Factor 
Emissivity Factor 



625 

89 (> 87) 

is-i±i:i 

FO II-III? (post-AGB) 
7750 
< 1-0 
135 ±5 
2.01 ±0.07 
^+0-9 

-0.8 



2.2 1 



B5±l V 

15,000 

4.0 

3.9±0.4 
5.9±0.8 



550 ± 50 

n o + O.l 
d -°-0.4 

0.475 
< 1 
£ 10 
0.3 
2.43 



HIPPARCOS (Pcrrvman et al. 1997) 
this work, Lissa uer*et"aL* 7T9*9'6*^ — ** 
this work 

this work 

this work, Castelli (1978 ) 
this work, Castelli (1978 ) 
this work 
this work 
this work 

this work 
Cox (2000) 
^ox (2000) 
"Cox" (2000) 
Cox (2000) 



Lissauer et al. 




Lissauer et al. 



3.2. The F Star 

We utilized the T = 7750 K, logs = 1-0 (|Castellil 
Il978f ). solar abundance model spectrum from the grid of 
iCastelli fc Kuruc3 (20041^1 to represent the F star com- 
ponent in e Aur, which has the appearance (more on that 
later) of an F0 supergiant star. 

We app lied limb- darkening to the model spectrum 
using the iVan Hammd (|1993l) relation, which is linear 
in form with wavelength-, temperature-, and gravity- 
dependent coefficients. The optical-near-IR region of the 
SED is dominated by this component, and we scaled the 
limb-darkened model spectrum to match the J -ban d pho- 
tometric point from iTaranova fc Shenavr in (200l[) (see 
Table [T]). This requires a radius of Rf — 135 Rq, re- 
sulting in an angular diameter of 2.01 mas. The radius 
can only be changed by ±5 Rq, corresponding to ±0.07 
mas in angular diameter, without significantly worsening 
the match to the J point. Our model angular diameter 
is somewhat smaller than t he va lue of 2.27 ± 0.11 mas 
measured by iStencel et al.l (|2008T ) in the K band, but is 
in agreement within 2a. 

The expected surface gravity of a "normal" F0 super- 
giant with radius of 135 Rq and mass of 15-20 Mq is 
log g ~ 1.5, which is one-half dex larger than the value se- 
lected for our model spectrum for this component. How- 
ever, as we will demonstrate below, the mass of the F 
star in e Aur is significantly lower than that of a nor- 
mal F supergiant. In fact, a more appropriate value for 
the gravity would be log g = 0.5; unfortunately, model 
spectrum calculations are apparently not typically made 
for stars that look like F supergiants but have masses an 
order of magnitude smaller, so log g = 1.0 is the lowest 
gravity model available. In any case, the gross differences 
in the model spectra between log g values of 1.0, 1.5, and 

12 The grid of model spectra is availa ble at 
http: //wwwuser . oat . ts . astro . it/castelli/grids .html and 
ftp : //ftp . stsci . edu/cdbs/gr id/ck04models/ 



2.0 can be compensated for by changing the stellar radius 
by only a few per cent, comparable to the determined ±5 
Rq uncertainty, so we expect no significant difference in 
the overall appearance of the model SED. 

3.3. The B Star 

Shortward of rs 0.15 /mi, the F star spectrum drops 
off sharply, and the B star spectrum dominates the SED 
of e Aur. This far-UV wavelength region (see Figure 
[2|) provides the primary constraint on the spectral type 
(and, hence, mass and radius) of the B star. Using the 
interstellar dereddening of the observed SED (see Sec- 
tion l3.1[) . a B8 V star (represented by the limb-darkened, 



T= 12,000 K, 



4.0, 



solar abu ndance model spec- 
initially provides a 



trum from [Castelli &: Kurucd 12004 
good match to the observed SED in this region, with 
the exception that the Ly-ct absorption line in the model 
spectrum is significantly broader and deeper than ob- 
served in the HST Goddard High Resolution Spectro- 
graph (GHRS) spectrum of e Aur. 

However, the B star is completely embedded in the 
dusty disk, which is viewed close to edge-on and has a 
significantly larger thickness than the radius of any B 
type main sequence star (see Section |3T4"|) . Consequently, 
we must consider that the observed spectral contribution 
of the B star will necessarily have been modified by pas- 
sage through the disk. There are two likely effects. The 
first is that the observed spectral contribution of the B 
star will be overall fainter compared to the unobscured 
case; this is essentially the same thing that happens to 
the F star during eclipse, but affects the B star at all 
times. In the example mentioned above, the B8 V star 
is only a good match to the observed SED if the disk is 
completely transparent, which, among other inconsisten- 
cies, would preclude the possibility of the disk eclipsing 
the F star. 

At optical-near-IR wavelengths, the eclipse of e Aur 
has been observed to be approximately gray (e.g.,[Kopal 
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Figure 2. Expanded view of the ultraviolet wavelength region of Figure [T] 



11971b iLissauer et al.1 1199 6'! . This is taken as evidence 
that the dust grains in the disk around the B star must 
be relatively large, a conclusion that is also supported 
by the lack of broad dust emission features in the Spitzer 
IRS spectrum of e Aur (see Figure EH), which implies 
dust grains of size > 10 mm (see iMorales et all 120091 : 
iNatta et all 120071 : iD'Alessio etaTl I2006L and references 
therein, for discussions of the effect of dust grain size on 
observed spectral features). We can simulate this effect 
by applying a simple scaling factor (which we will refer to 
as the disk "transmissivity" factor) to the B star model 
spectrum. A transmissivity factor of 1.0 would corre- 
spond to a completely transparent disk, such that the 
light from an object inside or behind the disk would be 
unaffected. On the other hand, a transmissivity factor of 
0.0 would correspond to a completely opaque disk, such 
that the light from an object inside or behind the disk 
would be completely blocked. For example, if the disk 
attenuates the spectral contribution of the B star by 50% 
(a transmissivity factor of 0.5), then the observed SED 
could still be matched by the presence of two B8 V stars, 
cor responding to the close bina ry B s tar model suggested 
by ILissauer fc Backmanl (|1984f ) and lEgglcton fc Pringlel 



19851) . 

The second likely effect on the observed B star spec- 
trum introduced by the circumstellar disk is an addi- 
tional reddening. The difficulty of distinguishing be- 
tween the observable effects of inter stellar and circum- 
stellar material is often no t ed (e.g... iHempel fc Schmittl 
l200l IHolweeer et all [19991 IWood et al.l 119971) . In the 
case of e Aur, we have the relative advantage that the 
interstellar reddening can be determined in the optical 
region of the spectrum. In this region, the F star domi- 



nates the light from the system, and is unaffected by the 
dusty disk (outside of eclipse). However, this still leaves 
the question of how to account for the circumstellar red- 
dening of the B star. 

To first order, additional circumstellar reddening of the 
B star might be blamed on the fact that the effect of red- 
dening is an order of magnitude greater in the UV than 
in the optical. Consequently, the "grayness" of the disk 
could begin to break down at the short wavelengths at 
which the B star spectrum dominates the SED. There 
is som e disagreement in the literature whether this ef- 
fect is (iParthasarathv fc: Fruehl 119861: lAltner et aLlll984D 
or is not (|Ferluga fc Hackiri985l ) supported by UV obser- 
vations with IUE. 

Another plausible scenario is that the hot B star is sur- 
rounded by a localized shell or torus of gas and small dust 
grains created by the sublimation and photo-spallation of 
the larg e dust grains that com prise the bulk of the disk 
(also see Cha pman et al.|[l983D . This material would pro- 
duce an additional non-gray attenuation of the B star's 
spectrum, without affecting the bulk properties of the 
entire disk. To account for this, we applied additional 
dereddening, corresponding to a larger value of E(B- 
V), to just the B star-dominated porti on of the SED. 
We used the iFitzpatrick fc M assa ( 20071) reddening law, 
thereby making the implicit assumption that extinction 
caused by the localized material around the B star has 
the same wavelength dependence as that of the interstel- 
lar medium (ISM). We note that this scenario is consis- 
tent with the suggested presence of a "hole" (probably 
a decrease in optical depth rather than an actual ab- 
sence of disk material) at t he center of t he e Aur disk, 
which has been proposed by Wilson ( 19711 ) and others to 
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Figure 3. Expanded view of the mid-infrared wavelength region of Figurc[T] 



explain the mid-eclipse brightening observed during the 
1982-1984 eclipse. We do not further consider the hole 
since we are not concerned with the eclipse behavior, ex- 
cept inasmuch as the typical eclipse depth implies that 
a specific fraction of the F star's light is blocked - see 
Section [3~6l 

The additional dcreddening correction was applied to 
the UV region of the SED as follows. For A < 0.15 
/Ltm, the SED is completely dominated by the B star 
component, so we simply applied a larger dereddening 
(i.e., interstellar + circumstellar) to these data than 
to the rest of the SED data (i.e., interstellar only). 
We applied a combined dereddening corresponding to 
E(B-V) = +0.45 (reddening steps of +0.05 in E(B-V) 
between E(B— V) — +0.40 — h-0.70 were considered). For 
A = 0.15-0.33 yitm, the SED is dominated by neither the 
B star nor the F star component. In order to construct 
the final dereddened data in this wavelength range (con- 
sisting of the WE spectra) , we first created two versions 
of the data dereddened using E(B — V) = +0.38 (in- 
terstellar only) and E(B - V) = +0.45 (interstellar + 
circumstellar), and summed them after weighting by the 
relative contributions of the stellar components to the 
total SED. For example, at A ~ 0.15 /im the two stars 
contribute almost equally, whereas at A ss 0.3 /Ltm the B 
star component contributes less than 1% of the total flux 
density. 

In the case of two B8 V stars (as discussed above), 
in the presence of even a small amount of additional 
circumstellar reddening, corresponding to E(B-V)=0.40 
(i.e., an increase of only 0.02 mag compared to the in- 
terstellar reddening alone), the disk transmissivity must 
be increased to 0.7 to preserve the match to the ob- 



served SED. For reddening corrections of E(B-V)>0.45, 
the B8 V model continuum slope no longer matches that 
of the observed SED, with the latter being flatter than 
the former at wavelengths < 0.15 //m. Producing the 
best fit for these larger circumstellar reddenings, which 
matches the model to the SED at 0.15 /im, and pro- 
gressively underestimates the observed SED flux density 
at shorter wavelengths, requires a disk transmissivity of 
> 0.95 for E(B-V)> 0.45, respectively. This is, again, 
ruled out by the constraint that the disk cannot be com- 
pletely transparent and eclipse the F star. 

However, a single B5V star (the next available ear- 
lier spectral type, represented by the limb-darkened, 
T = 15, 000 K, log q = 4.0, solar abundance model from 
iCastelli &: Kurucd l2004) , combined with a disk transmis- 
sivity of 0.3, provides a good match to the UV region of 
the SED. This includes a much improved match to the 
width and depth of the observed Ly-a absorption line 
(see Figure [2]), which is significantly broader and deeper 
in the B8 V model, and increasingly too narrow and shal- 
low for B3 V and earlier spectral type models (see below). 

Spectral types of B3 V (the next available model tem- 
plate) or earlier for the B star are not viable solutions 
(even after considering additional attenuation and/or 
reddening by the disk) because as the B star becomes 
hotter, the continuum slope in this region becomes in- 
compatible with the observed (dereddened) SED. Specif- 
ically, while the observed SED continuum slope rises to- 
ward longer wavelengths up to 0.15 /im (at which point 
the B star ceases to dominate the SED), the continuum 
of a B3 V star is almost flat and the continuum of a 
B0 V star has started to rise slightly toward shorter wave- 
lengths. This holds true even after testing by the applica- 
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tion of implausibly large additional dereddenings. Hence, 
the B star spectral type is fairly tightly constrained to 
B5 V, with an uncertainty of approximately one spectral 
type in either direction. 

3.4. The Dusty Disk 

At a wavelength of s» 3-4 /im, the observed SED of 
e Aur begins to show a deviation from the F star spec- 
trum. This manifests as an IR excess such that, at 100 
/zm, the observed SED flux density is w 3 times brighter 
than that of the F star alon43- 

We must draw a distinction between the Spitzer data, 
which were obtained at orbital phases of « 0.8-0.95, 
and the other mid-IR data (from the ground and MSX), 
which were obtained at an orbital phase of ~ 0.5. The 
latter data are systematically somewhat brighter than 
the IRAC data in the s» 3-5 fim region, and show some 
deviations at longer wavelengths, as well. Some of this 
behavior, particularly the stochastic deviations at the 
longer wavelengths, can likely be attributed to the dif- 
ficulties in calibrating ground-based mid-IR photome- 
try However, the systematic difference at the shorter 
IR wavelengths suggests an actual difference in the char- 
acteristics of the cool compon ent, either since befor e the 
1982-1984 eclipse (when the iBackman et all 119841 data 
were obtained) , or within an orbital cycle since after the 
last eclipse. The latter scenario is consistent with view- 
ing the hotter side of the disk that is most irradiated 
by the F star in the data from <j) 0.5, compared to 
the case in the Spitzer observations, when the opposite, 
cooler side of the disk is most visible. This effect has 
been noted in the pa st (e.g.. lTaranova fc ShenavrTnll2001l : 
iLissauer et al.l[l996f ). and the relevant IR data shown in 
Figure []] are consistent with a suggested blackbody tem- 
perature of 800-1000 K (having peak flux density at 5-6 
fim) for the hot side of the disk. Future improvements 
to the disk model would benefit from a well-sampled IR 
data set that spans the entire orbital cycle of e Aur; of 
course, the long orbital period makes this difficult to ob- 
tain in practice. 

In the meantime, we have concentrated on reproduc- 
ing the IR excess delineated by the Spitzer data, since 
these data extend to the longest IR wavelengths in our 
data set, and have a homogenous, well-tested flux cali- 
bration. This IR excess can be reproduced remarkably 
well with a simple single-temperature (550 K) blackbody 
function, which has peak flux density at sa 9 /im. A \ 2 
minimization test applied to the data longward of 1 /im 
(excluding the data longward of 3 /im that were obtained 
near orbital phase 0.5 - see above) gives a best range of 
500-600 K for the temperature. This is the range within 
which the change in x 2 is less than 10% compared to the 
X 2 minimum (which occurs at Tbb = 550 K). Coinciden- 
tally, it is also the range within which the change in x 2 
in temperature steps of 25 K is always less than 10%. 

13 We note that lAltenhoff eTall (1994ft measured a 250-GHz 
(1200- jUin) flux density for e Aur of 9±2 mjy, which is more than an 
order of magnitude brighter than the Rayleigh-Jeans extrapolation 
of the F star model spectrum at that wavelength, fp 1200 = 0-8 
mjy. In fairness, this is also a factor of 3 brighter than the sum 
of the flux densities of the F star and blackbody models discussed 
here at 1200 /im, which is 1.3 mjy. This possibly implies that there 
are non-thermal emission processes contributing at extremely long 
wavelengths. 



ILissauer et al.l ()1996ft calculated a detailed model for 
the disk in e Aur, but for our purposes, it is not nec- 
essary to use a model that is so complex. We acknowl- 
edge that our model deals with the disk in terms of av- 
eraged bulk properties (uniform cylindrical volume with 
a uniform mass distribution) rather than more "realis- 
tic" characteristics such as a specific radial density pro- 
file, scale height, and so on. The difficulty of specify- 
ing the detailed physics of such a disk, in the absence 
of sufficient constrai nts, is clear from the dis cussions in 
ILissauer et al.l (|1996[ ) and lCarroll et al.l(|1991ft . However, 
it is also clear that the SED is reproduced very well by 
a parametrically simple model. Thus, the task of recon- 
ciling this simple model with detailed dust disk physics 
awaits even more detailed and comprehensive future ob- 
servational constraints (e.g., via inte rferometric imaging 
of the disk; iKloppenborg et al"1l2010ft . 

In its simplest form, our single-temperature blackbody 
model for the disk is parameterized by a single scal- 
ing factor which (in conjunction with the distance to 
e Aur) is related to the projected emitting area of the 
disk (hence, its size). Used in this fashion, this would 
correspond to a completely opaque disk that emits only 
from its visible projected surface area. However, as dis- 
cussed in Section 13.31 we have reason to assume that 
the dusty disk in e Aur must be partially transmissive. 
As discussed in Section 13.61 the completely opaque case 
provides only a useful limiting case for the estimation 
of physical parameters of the disk. This has led us to 
express the blackbody component scaling factor as the 
product of two parameters. The first of these relates to 
the projected surface area of the disk and the second is a 
disk "emissivity" factor. The latter parameter scales the 
blackbody flux to account for the emission of dust grains 
inside the partially transmissive disk that are visible from 
outside the disk and/or the fact that, while a true circum- 
stellar disk likely contains a range of dust temperatures, 
we are utilizing only a single-temperature model. This 
parameter goes hand in hand with the transmissivity fac- 
tor that we also assigned to the disk (see Section I3.3[) , 
which specifics how much of the flux of an object behind 
or inside the disk will still be visible through the disk. 

Both the transmissivity and emissivity factors are 
physically meaningful, but purely parametric in their ap- 
plication to our model. The former is constrained some- 
what by the value required to match the far-UV SED 
with a model B star spectrum. The value of the lat- 
ter, however, is determined almost solely by minimizing 
the x 2 °f the IR (A > 1 fim) region of the model. The 
caveat to this is that in general the emissivity would be 
1.0 for a completely opaque disk in which the flux con- 
tribution from the disk is determined solely by the pro- 
jected surface area - hence, dimensions and inclination 
- of the disk. For a more transmissive disk, the value of 
the emissivity factor will be larger than 1.0. 

3.5. Stellar Component Masses 

If the F star in e Aur is a bona fide F0 supergiant 
(luminosity cl ass I) , then it must have a mass of Mp w 15 
Mq or more (Cox 2000). The well known mass func tion 
of e Aur, / = 3.12 M (e.g., ILissauer et all 119961 and 
references therein) , then requires that the mass of the B 
star component be 13.7 Mq or more. This, in turn, would 
require the B star to be of spectral type Bl V or earlier, 
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which is excluded by the observed far-UV SED, which 
constrains th e spectral type to B5±l V (with a mass of 
5.9 ±0.8 M^lCoxl 120001). Am ong other arguments (e.g., 
see lLambert fc Sawvenfl986T h a massive evolved B star 
can be excluded as a plausible scenario by considering 
the small relative contribution of the B star to the SED. 
Compared to the factor of > 10 increase in stellar r adius 
betw een B5-8 stars of luminosity classes V and I ()Coxl 
2000) , such a star would be quite overluminous compared 
to the observed SED. Although some additional mass in 
the B star component can be attributed to the dust disk, 
it seems unlikely that > 8 Mq could be accounted for in 
this manner (see Section 13.61 for confirmation of this). 

If the B star is a B5±l V star, and assuming for the 
moment that the mass of the dust disk is negligible in 
comparison (i.e., <C 1 M ), then the known mass func- 
tion requires that the mass of the F star is 2.2^Qg M . 
This conclusively points to the identification of the F 
star in e Aur as a low mass po st-AGB star rather than a 
normal high mass supergiant ([Eggleton fc Pringlelll985l 
Saito et all 119871 lTakeuti|[l986l and ILambert fc Sawverl 
19861 have also suggested low mass F star models for 
e Aur, although only the last of these has explicitly ex- 
pl ored a post-AGB ident ification; also see the review 
in IGuinan fc D ewarf 20Q3). Additional support for this 
conclusion is provided in Section @] 

3.6. Disk Size and Mass 

There is considerable evidence that both the disk and 
the orbital plane of e Aur must be viewed close to edge 
on. For example, the flatness of the eclipse profile re- 
quires that if the eclipsing body is an approximately cir- 
cular disk, then it must be seen in projection at a high 
inclination (i.e., close to edge on). A circular disk viewed 
in projection at a low inclination (i.e., close to face on) 
would produce a significant variation in the geometrically 
obscured area of the F star as the eclipse progresses. In 
fact, this effect rather strongly con strains the inclinatio n 
to be in the range i > 87° (also see lLissauer et aL| [l996). 

In the limiting case where the inclination is i = 90°, 
the obscuration of the F star during eclipse must be pro- 
duced solely by the disk rim. We parameterize this with 
a disk thickness, expressed as the height measured from 
the disk mid-plane, ffdisk (i- e -: one- half of the actual geo- 
metric thickness of the disk) . During eclipse, the optical- 
near-IR flux density is reduced by «50%. To first approx- 
imation, this would result if the disk completely obscured 
50% of the projected area of the F star, which could be 
accomplished with a uniform disk having a thickness of 
0.5 AU (i.e., ffdisk = 0.25 AU). However, in order to 
reproduce the observed mid-IR SED, under the assump- 
tion that the disk is completely opaque, would require a 
disk radius of « 20 AU. 

If the masses of the stellar components in e Aur are 
M-p ss 2.2 and Mb ~ 5.9, then we can calculate the or- 
bital separation of the stellar components to be a sa 18.1 
AU, which clearly rules out a disk with a radius of 20 
AU. This orbital separation, in turn, allows us to esti- 
mate t he outer radius of th e disk (e.g., using equation 
2 from iLissauer et al.lll996l ). to be R disk w 3.8 AU. In- 
cidentally, the Keplerian orbital velocity at the edge of 
such a disk would be 37 km s _1 (starting from the gen- 
eral form of Kepler's Third Law, r 3 = P 2 M* and re- 



arranging to yield u ro t = 27Ta/ M*/r for units of dis- 
tance in AU, time in yr, and mass in Mq). This is 
comparable to the range of « 30-40 km s _1 inferred 
for the e Aur disk via o bservations of radial velocities 
of H and metal li nes (e.g. . IFerluga fc Mangiacapralll991l : 
ISaito et al.lll987t ILambert fc Sawverl 1 1 9861 ) . 

A disk with thickness of 0.95 AU (i.e., ffdisk = 0.475 
AU) and radius of 3.8 AU, viewed at an inclination of 
89°, is sufficient to produce a geometric obscuration of 
the F star's disk of 72% (i.e., 28% of the F star is com- 
pletely unobscured by the disk during eclipse). A disk 
transmissivity of 0.3 (see Section |3~3|) allows 30% of the 
flux from the obscured portion of the F star (i.e., 22% 
of the total flux) to be visible through the disk, result- 
ing in the 50% value constrained by the eclipse depth. 
Finally, an emissivity factor of 2.43 is required to pro- 
duce the model shown in Figure [T] (i.e., the observed 
flux of the slightly transparent disk is a factor of 2.43 
larger than would be observed if only the opaque outer 
surface of the disk contributed). Additional solutions at 
other inclinations (always > 87°) are possible for other 
combinations of disk height, transmissivity, and emis- 
sivity (e.g., Figure Q] shows the model produced by the 
disk parameters listed in Table [TJ however, with rdisk 
fixed at 3.8 AU and transmissivity fixed at 0.3, identi- 
cal eclipsed fraction and model SED are achieved using 
i = 88°, ffdisk = 0.725 AU, and emissivity of 1.54, or 
using i = 9 0°, ffd s k = 0.375, and emissivity of 3.41). 
As noted in iHuand (|1973l 119741 ) , the uniqueness of disk 
models for e Aur is problematic. 

Up to this point, we have assumed that the mass of the 
dust disk is negligible compared to the stars in e Aur. Let 
us now test the validity of that assumption, starting with 
comparisons to other astrophysical disks. For example, 
the dust disks around young (T Tauri) stars are found to 
have a sharply peaked distribution of total (gas -I- dust) 
mass centered on 0.01 Mq, with the majorit y of objects 
in the range 0.001-0.1 M (|Hartmannlll998[ ). However, 
these disks also have characteristic sizes (outer radii) of 
50-100 AU or more, implying that the e Aur disk, if 
similar in structure, would likely be more than two orders 
of magnitude less massive (much less if we consider only 
t he mass of dust ). 

iMorales et al.l (|2009f ) used Spitzer observations of 52 A 
and late-B type main sequence stars to estimate mini- 
mum dust disk masses of up to 0.6 M moon (« 2 x 10~ 8 
Mq). The minimum disk masses for most of their sam- 
ple were several orders of magnitude smaller. They note 
that the total disk masses (i.e., including larger items 
that evade detection), assuming parent body planetesi- 
mal sizes of 1 km, could be a factor of 10 4 larger (i.e., 
~ 10~ 4 M©). They modelled these disks as annular 
rings around the parent star, and found a median ra- 
dius of 11.9 AU (range of 5-93 AU) with a radial spread 
of ~ 3-40 AU. This would present a surface area larger 
by a factor of at least ~ 5 than th e disk in e Aur , so we 
would expect the disk masses from [Morales c t al.l (|2009f ) 
to overestimate the e Aur disk mass, if they are otherwise 
comparable in structure. 

Another estimate of the e Aur disk mass can be made 
by comparison to the circumbinary dust disk that pro- 
duces the mid-IR excess in the cataclysmic variable 
V592 Cassiopeiae. This disk, which is itself the most 
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massive dust disk yet known for a cataclysmic variable, 
has bee n modeled to conta in 2.3 x 10 21 g (« 10 -12 Mq) 
of dust (jHoard et al.ll2009D . If we scale up the V592 Cas 
dust disk to the same volume as the e Aur disk, then the 
mass of dust in the latter would be ~ 6 x 10 -5 Mq. It is 
not clear if this scaling up exercise is appropriate, since 
the morphology and, hence, presumed mass distribution 
of the circumbinary disk in V592 Cas is quite different 
(geometrically very thin compared to its radius, a la the 
rings of Saturn) compared to the e Aur disk (which is 
quite thick relative to its radius, almost toroidal). How- 
ever, the binary star hosts of these disks are similar in 
all but scale, so we might expect a scaled comparison 
of their dust disks to still provide some useful limiting 
values. 

From opacity argum ents, Hinklc fc Simonl (|1987l ) and 
iBackman et al.l (|1984f ) estimated disk masses of ~ 10 ~ 6 
Mq (gas only) and ~ 10~ 7 Mq (gas and dust), respec- 
tively. We can attempt to further quantify an upper 
limit to the disk mass in e Aur in a similar fashion. 
We start by assuming that the dust in the disk is com- 
prised solely of spherical silicate grains with density of 
3.0 g cm -3 and radius of 10 fj,m. If we imagine looking 
through a 1 cm 2 column of such grains, then on aver- 
age [(1 cm)/(0.001-\/2cmgrain _1 )] 2 = 5x 10 5 dust grains 
per cm 2 would be required to produce complete geomet- 
ric obscuration along the length of the column (i.e., any 
photon traveling along the column will inevitably be in- 
tercepted by a dust grain). The factor of y2 provides for 
the round cross-sections of the spherical dust grains to 
overlap in order to completely fill a 1 cm 2 square column 
cross-section. If the column length is of the order of the 
disk radius, then the mean dust grain number density is 
8.8 x 10 -9 cm -3 , corresponding to a total mass of dust 
in the disk of « 8 x 10~ 9 Mq. Of course there could 
be more mass than this - we cannot discriminate past 
the point at which the disk would be completely opaque 
due to the column density of dust grains. However, the 
fact that we do see the B star inside the disk implies 
that (as already assumed in our model), the disk is not 
completely opaque, so this mass estimate is necessarily 
an upper limit. 

The mass of a spherical dust grain with a fixed density 
increases in proportion to r^ nst , while the cross-section 
of the grain increases as r 2 ust . Thus, increasing the grain 
radius by a factor of 10 decreases both the requisite col- 
umn density to achieve total obscuration and the total 
number of dust grains in the disk by a factor of 10 2 , but 
increases the total mass of the grains only by a factor 
of 10. In this manner, we can estimate that if the to- 
tal mass of dust in the disk is ~ 1 Mq (and the disk is 
completely opaque), then the "dust grains" would have 
to have implausibly large radii, rdust ~ 1 km. 

Therefore, based on all of the comparisons and calcu- 
lations described above, it appears safe to say that the 
mass of the dust disk in e Aur is <C 1 M Q , with the range 
of estimates topping out at around one-hundredth of a 
per cent of a solar mass. There is presumably also gas 
in the disk, but even at the "standa rd" ISM gas-to-dust 
mass ratio of -100:1 (jTielensI 12005). the mass contribu- 
tion from the gas will still leave the estimated total disk 
mass at well under 1 M . Plausibly large dust grains 
(e.g., centimeter-scale "pebbles") are expected to have 



distinct sub-mm emissivity properties (pe rhaps related to 
the 9 mJy 250-GHz flux density found bv lAltenhoff et all 
1994), and deserve further study at long wavelengths. 

4. EVOLUTIONARY MODELS FOR EPSILON AUR 

Just after the last eclipse event of eAur, IWeb bink 
(|1985|) presented a few model scenarios for the evolu- 
tion of the binary. Given the poorly known distance at 
the time and its MK spectral classification, e Aur was 
considered an F supergiant based on its spectral char- 
acteristics (e.g ., narrow lines) and effective temperature. 
ICastellil (|1978f ) performed a spectroscopic fine analysis to 
derive log g = 1 , indicative of a low surface gr avity con- 
sistent with a very large star. IWebbinkl (|1985l ) assumed 
three possible absolute luminosities for the supergiant, 
bracketing the full range of the distance estimate at the 
time. His models for e Aur initially included a pre-main 
sequence star, making the binary extremely young, but 
then he quickly argue s against such a scenario on a num- 
ber of grounds. Next. IWebbinkl (| 19851) considered a num- 
ber of shell and core burning scenarios, all more or less 
suggestive that e Aur was a post-main sequence star lo- 
cated in one or another "loop" of its evolution crossing 
the H-R diagram. Some of these loops involve significant 
mass loss from the star during an earlier stage. 

While the general idea of the presence of a post-main 
se quence evolved F star in e Aur is not a new idea (e.g., 
see lLambert fc Sawverll986i ) , we have proposed here that 
the apparent F supergiant star is, in fact, a bright post- 
asymptotic giant branch (post-AGB) star that started 
on the main sequence as a > 7 Mq star. The ob- 
served 12 C/ 13 C ratio of ~ 10 from CO line observa- 
tions (Hinklc fc Simonlll987|) and possibly e levated bar- 
ium abundance of 2 x solar (jCastellil 119781 ) are indica- 
tive of AGB thermal dredge-up and s-process enhance- 
ment appropriate for a post-AGB star. W e note that the 
elevated barium abundance result from (Castclli 1978) 
would benefit from confirmation using updated oscillator 
strengths, and in comparison to line strengths of other 
s-process elements. We believe that the mass of the F 
star is currently near 2.2 Mq, while its size is that of 
a supergiant (135 Rq). The constraints on the mass 
and radius of the F star, which include comparison with 
multi-wavelength ( i.e., the SEP presen ted here) and in- 
terferometric fe.g.. iStencel et aLll2008D observations and 
the HIPPARCOS parallax distance, as well as known 
kinematic properties of e Aur, are discussed in detail 
above. We justify the plausibility of the proposed low 
mass, post-AGB F star below. 

The study of post-AGB stars is a developing field in 
stellar evolution researc h, as more examples are identi- 
fied observationally (see ISzczerba et al.l 120071) . Among 
the few single star evolutionary models that come 
close to the post - AGB phase, are those presented by 
iMarigo fe Girardil (|2007l ). which include thermal pulse 
tracks for stars up to 5 Mq and a range of metallic- 
ity. The inherent instabilities limit the full exploration 
of these extreme late phases of evolution, but such work 
provides guidance for interpreting the F star in e Aur. 
Given the observational constraints on the F star, namely 
T eff , R/Rq, logi = 4.7, Z ~ Zq and the - 100 d 
q uasi-period of photospher ic variations, the 5 Mq tracks 
in IMarigo fc Girardil (|2007l ) come closest to approaching 
these values and have evolutionary timescales approach- 
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ing 10 5 yr. By extension, we could place the progenitor 
of the F star in the 6-8 Mq range, making it a cand i- 
date for so-called super- AGB status (|Smartt et al.ll2009l) . 
Tidal interaction between components in e Aur might ac- 
count for the inflated size of the F star and some of the 
out-of-eclipse light variations. 

The effective temperature and luminosity of the F star 
component in e Aur matches well with an initially 7 
Mq, post- AGB star early in its evolution away from the 
AGB, probably after being "born-again" by a thermal 
pulse due to rapid core burning. However, the evolution 
tracks in the H-R diagram for the initial movement away 
from the AGB and following a thermal pulse are indis- 
tinguishable, so we can not tell where in its post- AGB 
evolution the F star is actually located. In any case, 
these "loops" on the H-R diagram are extremely rapid 
for su ch a star, with durations of only a few thousand 
years (jBloeckerl 119951 ). Ancient star catalogs, however, 
do not imply much change in the appearan ce of e Aur on 
the t i mescale of recorded hum an history ()Carroll et al.l 
[1991: IGuinan fc Dewarfl 120021) . 

The known mass function of e Aur requires that the F 
star currently has a mass of ~l-3 Mq for a ~5-7 Mq 
B star (where the range in B star mass represents the 
constrained range of its spectral type, B5±l V). By im- 
plication, during the past few 100,000 years, the F star 
must have lost ~ 5 or more solar masses of material dur- 
ing its AGB and post-AGB evolution. This high mass 
loss amount argues for the star being in its first or even 
second thermal pulse loop tour of the upper H-R dia- 
gram. The F star was the initially more massive star in 
the e Aur binary system, starting on the main sequence 
as an early B star. This assignment gives e Aur an age 
of approximately 60-80 million years since the ZAMS. 

The additional complic ation of binary star evolution 
requires further work (see |Eldridgc et al. 2008]), but as- 
sessing the F star to be in a transitional state (blue 
loop) in its evolution makes an attractive and testable 
hypothesis. The implied 10 4 yr or longer timescales 
could provide ample time to transfer and lose mass in 
and around the binary, but not allow the resultant disk 
to dissipate, nor the photosphere to chemically modify 
due to hot bottom burning. The B5 V star and its 
surrounding dust disk are likely capture sites for some 
of the mass lost by the present-day F star, but not as 
much as one might think. As the precursor F star ex- 
panded after its main sequence life was over, the binary 
nature of the system would have allowed matter to flow 
through not only the inner Lagrange point (LI, located 
between the two stars) but also a greater flow would 
have escaped through the outer Lagrange point (L2), as 
model ed and discussed in lEggleton fc Kiseleva-Eggletonl 
(2002) in relation to Algol-like binaries. The issue of 
large amounts of mass loss in close binaries, especially 
the idea of non-conservative mass loss through the outer 
Lagrange points, is famously known as the Algol Paradox 
and was treated in detail in iKopall (|1978| ). 

If the mass lost from e Aur leaked out of the L2 side 
of the binary, then it might remain "hidden" in a thin, 
flat ring of measurable extent. Observers might look for 
evidence of faint, tell-tale circumbinary emission from 
such residual material around e Aur. Taking the ad- 
vanced post-AGB evolution nature of the primary star, 



we believe that the "supergiant" size is due to contin- 
ual mass loss via a slow expanding photosphere heading 
toward shell ejection as in a planetary nebula or simi- 
lar type of object (e.g.., Sak urai's object or FG Sagittae; 
lLawlor fc MacDonald 120 3). The inferred F star radius 
of 135 Rq is consistent with these stars and the models 
produced for them. 

5. CONCLUSIONS 

We have analyzed an unprecedented SED of e Aur, as- 
sembled from observational data spanning three orders 
of magnitude in wavelength. In conjunction with con- 
straints provided by other published information about 
this enigmatic binary star, such as its orbital dynamics, 
we can conclude that: (1) the F star component is a low 
mass post-AGB object, (2) the B star component is a 
B5±l V star, and (3) the dusty disk is a partially trans- 
parent, low mass disk of predominantly 10 fim or larger 
grains. The identification of the F star in e Aur as a 
normal high mass F supergiant is simply no longer ten- 
able as a plausible scenario. The requisite mass of the B 
star plus disk cannot be made to satisfy the well-known 
mass function for this system (without invoking exotic 
scenarios involving compact, non-luminous sources of 5- 
10 Mq of additional mass that do not produce X-rays 
- see Section [1]), while at the same time satisfying the 
constraints on the luminosity and SED by the data we 
have assembled here. From an evolutionary standpoint, 
the F star must have been initially the more massive of 
the stellar components in the binary, and some or all of 
the mass in the dusty disk around the B star may have 
been transferred from the F star precursor. The bulk of 
the mass is likely to have escaped from the binary during 
the evolution of the F star precursor, and might be vis- 
ible in sensitive observations at far-IR wavelengths. As 
a post-AGB object, the F star is in a relatively rapid 
transitional phase of stellar evolution, and we should ex- 
pect significant changes in the appearance of e Aur after 
the next few thousand to tens of thousands of years. In 
the meantime, observational studies of the out-of-eclipse 
disk at wavelengths greater than 30 (im and via interfer- 
ometric imaging are to be encouraged. 
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